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Abstract 

We recalculate the two loop corrections in the background heat bath using 
real time formalism. The procedure of the integrations of loop momenta with 
dependence on finite temperature before the momenta without it, has been 
followed. We determine the mass and wavefunction renormalization constants 
in the low temperature limit of QED, for the first time with this preferred order 
of integrations. The correction to electron mass and spinors in this limit is 
important in the early universe at the time of primordial nucleosynthesis as well 
as in astrophysics. 
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1 Motivations 

The quantum field theory at finite temperature, Thermal Field Theory (TFT) 
has been used for more than 50 years. The techniques for calculations in many 
body systems were initially developed in condensed matter [1 ] which are now 
also used to describe a large ensemble of multi-interacting particles in a ther- 
mal background. Since the late seventies, TFT was used to study the phase 
transitions in cosmology and quantum field theories. It was adopted in particle 
physics to study the many particle systems. These methods are used not only 
in particle physics but also study the various aspects of nuclear matter and 
plasma physics. The TFT methods are extensively used to describe the phase 
transitions due to symmetry breaking after the Hot Big Bang and in tracing the 
history of the early universe. 

The main applications of interest in this context are in: 

Cosmology: The early universe provides a very good example in the studies 
of hot plasmas. When the Big Bang occurred, the Universe was impossibly hot 
and dense. It rapidly expanded and cooled. At t = 200 sec, T = 10 8 -?T, it 
was cool enough for neutrons and protons to combine to form Deuterium, then 
Helium and traces of Lithium (primordial nucleosynthesis). For the next few 10 5 
years it was too hot for electrons to form atoms. The universe was filled with 
hot plasma of electrons and nuclei, bathed in photons constantly interacting 
with both, like the interior of a star. 
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Astrophysics: There is a series of different types of fusion reactions in stars 
leading to luminous supergiants. When helium fusion ceases in the core, gravi- 
tational compression increases the core's temperature above 6 x 10 8 -£T at which 
carbon can fuse into neon and magnesium. As the core reaches 1.5 x 10 8 JsT, oxy- 
gen begins fusing into silicon, phosphorous, sulfur, and others. At 2.7 x 10 8 K, 
silicon begins fusing into iron. This process essentially stops with the creation 
of iron and a catastrophic implosion of the entire star initiates. When the high 
mass stars exhaust their He fuel they have enough gravitational energy to heat 
up to 6 x 10 8 K. Cores of neutron stars, red giants and white dwarfs are composed 
of extremely dense plasmas (p = 10 6 -10 15 g/cm 3 , ). The neutrinos and axions 
emission rates in these stars require TFT [2] . A tremendous amount of energy 
is released in a supernova. The only supernova in modern time, visible to the 
naked eye, was detected on Feb. 23, 1987 and is known as S7V1987AIt emitted 
more than 10 10 times as much visible light as the Sun for over one month and 
temperatures as high as 2 x lO^K were reached. 

Heavy Ion Collisions: The quark gluon plasma is the form of matter at 
transition temperatures T c = 100 — 200 MeV. The hot and dense environment 
in quark gluon plasma and the studies of its prospective reproduction in nucleus- 
nucleus collisions require the TFT methods for detailed explanation. With the 
increased feasibility of creation of quark gluon plasma in heavy ion collisions, 
the methods developed in this theory got their specific relevance in QCD at 
finite temperature as well. 

We are specifically interested in the first two applications here. 

2 Finite temperature effects 

The main idea of TFT is to use the approach of the usual quantum field theory. 
Matsubara [3] was the first who developed thermal field theory by incorporating 
a purely imaginary time variable in the evolution operator. In Euclidean space 
the covariance breaks and time is included as an imaginary parameter. The 
imaginary time domain is finite and periodic because of which the energy inte- 
grations are converted into summations over the discrete Matsubara frequencies. 
The presence of discrete energies along with the particle distribution functions 
destroys the covariance of the theory. 

The important contributions by Schwinger [4] , Mills [5] , and Keldysh [6] led 
to the development of a formalism based upon the choice of a contour in the 
complex plane. This is called the real time formalism. In the real-time for- 
malism, an analytical continuation of the energies along with Wick's rotation 
restores covariance in Minkowski space at the expense of Lorentz invariance. 
The breaking of Lorentz invariance leads to the non-commutative nature of the 
gauge theories [7] . The covariance is incorporated through the 4-component ve- 
locity of the background heat bath = (1, 0, 0, 0).ln a heat bath the particles 
are in constant interaction with the thermal surroundings. Implementing these 
interactions is straightforward as is done in vacuum field theory. The tempera- 
ture is included through the statistical distribution functions of the particles. 
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Umezawa and coworkers [8] independently worked on a different approach 
called Thermo- Field Dynamics that also gives the same results. In this formal- 
ism, the propagators are taken in the form of 2 x 2 matrices. Field theory at 
finite temperature is renormalizable, if the vacuum theory is so since the pres- 
ence of the Boltzmann factor in the thermal corrections cuts off any ultraviolet 
divergence. Choosing the suitable counter terms as in vacuum can eliminate 
them. The infrared divergences are inherent in almost all perturbation theories, 
whether at zero or finite temperature. KLN theorem [9] demonstrates that sin- 
gularities appearing at intermediate stages of the calculation cancel in the final 
state physical result. 

Quantum Electro Dynamics (QED) is the simplest and most successful gauge 
theory. The behavior of QED at finite temperatures serves as a model for the 
determination of background effects in other physical theories - the elcctroweak 
theories as well as Quantum Chromo Dynamics. In the real time formalism, the 
tree level fermion propagator in Feynman gauge in momentum space is [10] 

Sfi (P) = 0/ - m) [ 2 _ * - 2ir5(p 2 -m 2 )n F (E p )] , (1) 

P 771/ ~\ IS 

where 

M^ P ) = e (P i) + 1 . (2) 

is the Fermi-Dirac distribution function with /3=^. The boson propagator is 
taken as 

D^( P ) = [ w ^--2n6(k 2 )n B (k)}, (3) 

with 

= eW ,!)_ r (4) 



3 One loop corrections 

At the one loop level, Feynman diagrams are calculated in the usual way by 
substituting these propagators in place of the usual ones in vacuum. The Lorentz 
invariance breaking and conserving terms remain separate at the one loop level 
since the propagators comprise temperature dependent (hot) terms added to 
temperature independent (cold) terms. This effect has been studied in detail 
and established at the one- loop level [11]. The renormalization of QED in this 
scheme [12] has already been checked in detail at the one loop level for all 
temperatures and chemical potentials. 
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The thermal background effects are incorporated through the radiative cor- 
rections. In finite temperature electrodynamics electric fields are screened due 
to the interaction of the photon with the thermal background of charged parti- 
cles. The physical processes take place in a heat bath comprising hot particles 
and antiparticles instead of vacuum. The exact state of all these background 
particles is unknown since they continually fluctuate between different configu- 
rations. The net statistical effects of the background fermions and bosons enter 
the theory through the fermion and boson distributions respectively. 

The electric permittivity and the magnetic susceptibility of the medium are 
modified by incorporating the thermal background effects. At low temperatures, 
i.e., T « m e (m e is the electron mass), the hot fermions contribution in back- 
ground is suppressed and only the hot photons contribute from the background 
heat bath. The vacuum polarization tensor in order a does not acquire any hot 
corrections from the photons in the heat bath. This is because of the absence 
of self-interaction of photons in QED. 

The thermal mass is generated radiatively. The mass shift that enters physi- 
cal quantities acts as a kinematical cut-off, in the production rate of light weakly 
coupled particles from the heat bath. The effective mass corresponds to the fact 
that in the heat bath, the propagation of particles is influenced by their contin- 
uous interactions with the medium. 



4 Higher order corrections 

The higher order loop corrections are required to get predictions on perturba- 
tive behavior at finite temperature. At the higher-loop level, the loop integrals 
involve a combination of cold and hot terms which appear due to the over- 
lapping propagator terms in the matrix element. In such situations, specific 
techniques are needed, even at the two loop level, to solve them. Higher loops 
get analytically even more complicated. In the hot terms there appear overlap- 
ping divergent terms. The removal of such divergences is already shown at the 
two- loop level [13] for electron self energy. 

We restrict ourselves to the low-temperatures to prove the renormalizability 
of QED at the two-loop level through the order by order cancellation of sin- 
gularities. The results have been shown [14] to depend on the order of doing 
the hot and cold integrations. The justification of this specific order is the fact 
that the temperature dependent part corresponds to the contribution of real 
background particles on mass-shell and incorporates thermal equilibrium. The 
breaking of Lorentz invariance changes these conditions for the cold integrals. 
We have checked that the renormalization can only be proven with the preferred 
order of integrations, i.e., if covariant hot integrals are evaluated before the cold 
ones. At the higher loop level the vacuum polarization contribution is non zero, 
even at low temperature. The calculations are simplified if the temperature de- 
pendent integrations are performed before the temperature independent ones. 
The temperature independent loops can then be integrated using the standard 
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techniques of Feynman parametrization and dimensional regularization as in 
vacuum [15]. 

The second order in a corrections to the electron self energy at low temper- 
ature has been calculated to be equal to 



S /3 (p) = E^(p,T = 0)- 5 -3 
where the singularity 



- £ [{jf+ 6m)I A - 21] - 3(j/ + Am)! A - 4 7 ^ - 



8irT 2 
3m 2 
(5) 



G/-m) 



1 1 47T/X 2 

- = - - 7-M — 2"), 

is the usual ultraviolet divergence in vacuum field theory in MS bar scheme 
of renormalization with 7 as the Euler- Mascheroni constant and 



f dk . 
I A = 8tt / —n B {k), 



(6) 



7° 



2tt 3 T 2 , 1 + t; 
o — m . 

3E 2 v l-v 



(7) 



I.p 



3t2 



2tt j T 



fin 



1 + v 
l-v 



2t>], 



(8) 



with t> = — . The divergences of the form I a and - in Eqs. (6) and (7) cancel 
on addition of the appropriate counter terms. The finite contribution then leads 
to the physical mass of electron: 



m phys — m 



2anT 2 Aa 2 T 2 



3m 2 



3to 2 



(9) 



5 Results 

The electron mass and wave function renormalization can be obtained from the 
two loop self-energy for electrons calculated in the previous section. From Eq. 
(10) the change in the electron mass due to finite temperature up to the order 
a 2 relative to the cold electron mass is 



5m 
m 



T 2 



(air + 4a 2 ) . 



(10) 
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Figure 1: 
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Fig. 1 



The renormalizability of the self mass of electron is rechecked through the 
order by order cancellation of singularities at both loop levels. It can also be 
noted that the second order term is much smaller than the first order term. 

Using the standard procedure, the wave function renormalization constant 
comes out to be 

In Fig. 1, a plot of Z^ 1 vs T/m is given for 0(a 2 ) to demonstrate the 
affect of temperature in two-loops. The results here and in Rcf. [14] are an 
explicit proof of the renormalizability of QED up to the two-loop level. They 
also estimate the temperature dependent 

modification in the electromagnetic properties of a medium. This helps to 
evaluate the decay rates and the scattering crossections of particles in such a 
media. These results can be applied to check the abundances of light elements 
in primordial nucleosynthesis, baryogenesis and leptogenesis. If the background 
magnetic fields are also incorporated, then one can look for applications to 
neutron stars, supernovae, red giants, and white dwarfs. 

Figure caption 

Fig. 1 A plot of wavefunction renormalization constant Z^ 1 vs T/m at 
two-loop level. 
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